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ABSTRACT 



The sputtering of Cu fee (110), (100) and (111) surfaces by 1- 
to 10-keV A^ ions at normal incidence has been investigated by high- 
speed digital computer techniques. The interatomic repulsions are 
described by Bom-Mayer potentials. The shape of the sputtering 
ratio curves are in close agreement with the data of Magnuson and 
Carlston. Ejection patterns in accord with experimental observations 
are predicted. There is evidence that the sputtering mechanism is 
not the result of momentum focusing. 
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1. Introduction 



The process of ionic erosion of a solid surface is known as 
sputtering. As pointed out by Cheney, Rogers and Pitkin w and 
Harrison H , various theoretical models of the sputtering process 
have been developed since the 1850 T s when metal from the electrodes 
in gas discharge tubes was first observed to be transferred to the 
glass walls of the tubes. In 1923, Kingdon and Langpiuir [V] 
suggested a momentum transfer model, later elaborated upon by Lamar 
and Compton M • This theory was subsequently replaced by a 
thermal spike model, wherein a local hot spot caused evaporation 
(5,6^]. Later Keywell £ 7 ]] suggested a model based upon a succession 
of random collisions, while Harrison ^8-ljQ has proposed several 
models, including a statistical approach based on neutron diffusion 
theory, and Monte Carlo techniques. Pease O] developed a model 
based upon expressions from radiation damage theory. 

However, these early models neglected the ordered arrangement 
of atoms within the crystal lattice. Current theories are based 
upon a transfer of energy through collision processes, as developed 
by Silsbee [u] and Liebfried [i 6 ] . In addition to conducting a 
variety of experiments, Harrison [l7,18] combined these energy chain 
concepts with previous statistical theories into a radiation damage- 
statistical model [2,19] . 

Recently, workers in the field of radiation damage have employed 
digital computers to make theoretical studies. Robinson and Oen |2oj 
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have assumed that the '•knock-on 11 atom interacts through a series of 
binary (two-body) collisions, while Gibson, Goland, Milgram and 
Vineyard L*q. Erginsoy, Vineyard and Englert [22] , and Gay and 
Harrison £23] have used iteration techniques to numerically integrate 
the classical equations of motion of a many-body set of atoms. 

The present work uses the latter technique to investigate the 
sputtering of thi'ee orientations of a copper crystal subjected to 
bombardment by an argon ion. This particular combination of elements 
was chosen for two reasons; first, there is experimental data avail- 
able for comparison [24,25] and secondly, parameters to be used in 
the Born-Mayer potentials have been calculated [21,26,27] . 

Distinction should be made between "sputtering 11 and "transmission 
sputtering." As used in this paper, the former pertains to atoms 
ejected from the bombarded, or front, surface of a solid, while the 
latter refers to atoms leaving the back surface of a thin crystal. 
Although the possibility exists that an ejected atom will have one 
or more electrons stripped away in the sputtering process, this 
study makes no distinction between neutral atoms and ions. 

The sputtering ratio is defined as the number of atoms ejected 
per incident ion. Several aspects of sputtering have been studied 
experimentally [24,25,28,29,30, including the sputtering ratio as 
a function of primary (incoming ion) energy, incidence angle, crystal 
orientation and target temperature; the angular distribution of ejected 
atoms, and the energies of these atoms. This study is concerned 
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quantitatively only with the sputtering ratio as a function of 
monocrystalline orientation for primary energies in the 1- to 10-keV 
range and at normal incidence, and qualitatively with a description 
of the process whereby the crystal atoms are ejected from the lattice, 
and their angular distribution. 
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2 . The Model 



The model employed in this computer simulation program is the 
result of the combination of many ideas from previous work in the 
field of radiation damage. An early attempt by Oen, Holmes and 
Robinson [31] assumed binary collisions but excluded crystal lattice 
effects. This work suggested that neither the exponentially screened 
Coulomb potential nor the Bohr hard sphere potential were sufficiently 
accurate. A later model used by Robinson and Oen [ 20 ] placed copper 
atoms on their proper sites in a face centered cubic lattice. 

However, this model also considered collisions as binary and assumed 
that the target atoms did not interact amongst themselves. The Bohr 
potential and Born-Mayer potential with Gibson 2 parameters [2l] 
were compared. Gibson et<, al. [21] also studied effects in the fee 
copper crystal, but considered collisions as n-body while comparing 
several sets of the two parameters of the Born-Mayer potential. 

Similar work for bcc iron was done by Erginsoy et a al. £ 22 ]. More 
recently. Gay and Harrison employed a somewhat different n-body 

model. This project uses the Gay-Harrison model but with further 
modifications, vizs 

(1) The stationary immoveable atoms surrounding the crystallite 
have been removed. The object of this study is the ejection of atoms 
from the surface, and also, the lattice is stable before bombardment 
without these stationary atoms. 

(2) The crystallite size has been increased to contain the 
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sputtering event. The (110) sample contains 115 atoms, the (100) 
has 158, and the (111) contains 150. The (110) sample corresponds 
to ten planes in the x-direction, five planes in the 7 -direction, 
and nine planes in the z-direction. The (100) sample is 7 x 5 x 9, 
while the (111) is 9 x 5 x 17. The total external dimensions of 
these three crystallites are approximately the same (see appendix I). 
As work progressed it was readily determined which atoms were 
sputtered from the crystallite. Thus, by judicious choice of impact 
points, the crystallite may be decreased in size to reduce computer 
running time. This will be done in future studies. 

(3) A running report of the motion of every atom versus a 
final status report has been incorporated. 

(4) The impact area for the (111) sample has been chosen to be 
representative of a "repeat” crystallographic distance rather than 
just of the first two planes. The repeat distance is the number of 
layers that must be traversed to repeat the crystal configuration 
[ 24 ]. Thus the projections of the first two layers of the (110) 
sample, first layer of the (100) sample, and first three layers of 
the (111) sample are considered. Then symmetry shows that the 
rectangle for the (110), and triangles for the (100) and (111) 
orientations, will completely describe the crystallites if repeated 
(see figs. 1-3). Thirty-six impact points were assigned to the (110) 
representative rectangle, while 21 points were assigned to the (100) 
and (111) triangles (see fig. 4). 
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(5) There is no ’’target atom” as defined by Gay and Harrison, 
but the impact points are all near a "central" atom of the crystal- 
lite front surface. 

(6) The "time step multiplier" is 0.10 for all impact points 
except those in which a "head-on" collision with a front surface 
lattice atom occurs. In the latter case, a value of 0.05 was found 
to be necessary to minimize total energy error. 

In addition to n-body collisions, major ideas used by Gay and 
Harrison and retained in this program include the double iteration 
technique, lack of thermalizing motion, concepts of representative 
areas and time step intervals mentioned previously, the "lattice 
unit," disregard of interactions between the primary and lattice 
electrons (although provision is made for inclusion in future studies), 
use of Gibson 2 parameters in the Born-Mayer potential for copper- 
copper interactions, and omission of cohesive effects and restoring 
forces. The last point is justified by time considerations; the 
sputtering event is complete in a relatively short period, one-tenth 
to one-hundreth of the time necessary for these effects to commence 
[ 23 ] . 

Harrison, Carlston and Magnuson |j2€Q state that the repulsive 
exponential Born-Mayer potential function of the form V(r) = exp(A+Br), 
where B is always negative, appears to be adequate for most high 
energy collision dynamics in a lattice. Initially, the program was 
run using the Gibson 2 parameters in the Born-Mayer potential for 
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both the argon-copper and copper-copper interactions. These para- 
meters are not valid for the argon-copper system. The use of these 
constants resulted in a sputtering ratio versus primary energy curve 
of improper shape; the curve decreased with increasing energy with 
no sign of the experimental maximum around 3- to 5-keV. 

Subsequently, an argon-copper potential was employed for the argon- 
copper interactions, while the Gibson 2 was retained for the copper- 
copper system. This A-Cu potential is also of the Born-Mayer form 
and uses serai-empirical parameters calculated by Harrison et, al. 

from experimental data collected by Magnuson and Carlston [ 33 ] . 
The sputtering ratio curve obtained with this two potential approach 
exhibited a peak in the proper region, but the sputtering ratio 
magnitude was too low at high energies. To increase this ratio the 
A-Cu potential was '‘hardened** by increasing '’A*’ in the V(r) expression 
from the calculated value of 11.435 to 12.25. Figures 5,6 and 7 
result when the tabulated parameters are used in the Born-Mayer 
potential functions 



System 


A 


B (lattice unit) 


A-Cu 


12.25 


-7.597 


Cu-Cu 


10.024 


-9.197 



The computer program is a mixture of Fortran 60 and Fortran 
Symbolic language (see appendix II). The runs were conducted on a 
Control Data Corporation 1604 computer. Running time for each impact 
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point was five minutes for the (110) sample, nine minutes for the 
(111) and from ten to 15 minutes for the (100). Very generally, the 
program constructs the lattice, bombards the crystallite with a 
single primary, performs the iteration procedure to determine forces 
and displacements at the proper time steps, and stops when the 
total potential energy of all the atoms is less than 0.5 eV, or at 
a time step chosen by experience so that events of interest have 
already occurred. The output consists of the position and velocity 
of each atom with a total energy greater than thermal, assumed to 
be 0.25 eV (see appendix IV). Output is received when the run is 
25%, 50/S, 75% and 100% complete. 

At the end of the runs the totals of kinetic plus potential 
energies differed from the primary initial energies (the lattice 
has no energy until disturbed) by an average of 1%. This error is 
considered well within the uncertainties inherent in the model, 
especially the kinematic simplifications of the iteration procedure. 
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3. Results 



A. GENERAL DESCRIPTION OF THE SPUTTERING PROCESS 



When energy from the primary is transferred to lattice attorns 
through the n-body collision process, these atoms (secondaries) are 
displaced from their sites. The secondaries also interact in a 
series of collisions, and eventually some will gain momentum back 
towards the crystal surface. If any atoms have energies great 
enough to overcome the surface binding energy they will emerge from 
the crystal. Very generally, this is the sputtering process. 

Cheney et 6 al. £VJ have stated that sputtering is expected to 
be great when there is a large transfer of energy from the primary 
to the crystal lattice occurring reasonably near the surface, a 
rapid reversal of the momentum vector of the primary, a mechanism 
for the transfer of energy back to the surface, an easy escape of 
surface atoms and a comparatively large atomic volume or mass. 
Relatively open crystal structures permitting deep penetration of 
the primary results in a low sputtering yield. Thus, the yield is 
dependent upon the orientation of the crystal [ 24 ]. However, a 
crystal which has a history of radiation damage may not permit 
channeling |20,32j in open crystallographic directions, in which 



have suggested that the ejection of atoms along preferred directions 
is the result of momentum focusing chains. In light of this 
research some of these concepts must be modified. 



case the sputtering yield may actually increase. Other works 
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B. THE INDIVIDUAL ORIENTATIONS 



1. The (110) Sample 

Inspection of figure 8 reveals that most of the ejected atoms 
come from the region near the impact rectangle. A detailed analysis 
of the mechanism whereby a few of these atoms are sputtered will 
lead to an understanding of the sputtering process. In the following 
discussion '•low 11 energy refers to 1 keV and "high” energy refers to 
10 keV primary energy. 

Atom 14 is ejected most often. At low energy the impact points 
which cause sputtering of this atom are those in the upper right 
quarter of the representative rectangle. At high energy the impact 
points in the center of the rectangle cause it to sputter. The 
number of impact points which result in ejection is essentially 
independent of the primary energy. 

For all cases in which this atom is sputtered, it is driven 
towards the crystal interior by interaction with the primary. •‘Down*’ 
and '’left" motion is also imparted. The atom is repulsed by atoms 
lying in the second plane, no. 33 usually, and is reflected towards 
the front surface. This process imparts momentum to no. 33 towards 
the back of the crystal, which accounts for the relativly few times 
that it is sputtered. At 10 keV no. 14 gains enough momentum to 
interact also with atom no. 8 after rebounding from no. 33. 

Atom 15 is a major contributor to the 3 keV peak of the sputter- 
ing ratio curve. Impact points that cause this atom to sputter lie 
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in the right center portion of the representative rectangle at low 
energy, and the bottom center at high energy. No. 15 is ejected 
less often with increasing primary energy. Channeling, resulting 
in deeper ion penetration before interaction, occurs with greater 
primary energy, which explains this behavior. The mechanism which 
ejects 15 is dependent upon interaction at the front surface. The 
atom is initially directed to the right by the primary. The ensuing 
collision with no. 16 rebounds no. 15 out of the crystal. Addition- 
ally, at higher energies, the primary rebounds from no. 14 and 
ejects 15 directly by getting behind it before interacting strongly. 
Primaries entering at impact points in the left half of the rectangle 
do not interact strongly enough with no. 15 to impart the velocities 
necessary to collide with no. 16. 

Atom 37 is also a contributor to the sputtering ratio peak. 

The impact points responsible are in the upper left quarter of the 
rectangle, where the primary ’’bypasses” no. 14 and interacts mainly 
with second layer atoms, nos. 38 and 39 in particular. No. 38 is 
driven in and to the left in a direction to be behind no. 37 before 
approaching closely enough to ineract. Thus 37 acquires motion out 
of the crystal and towards the left. 

Comparable to atom no. 15, the impact points that sputter no. 

37 lie in a region where the primary is channeled at higher energies, 
with resultant less ion-atom interaction at the crystal front planes. 
No. 37, then, is also sputtered less at higher ion energies. 
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At low energy atom 19 is initially given a velocity up, in and 



towards the left. Interactions with no. 24 at the front surface 
and no. 43 in the second plane result in reflection towards the front. 
The impact region responsible for sputtering no. 19 is also the 
upper left quarter, with a few points around the perimeter of this 
area. With increased primary energy impact points where channeling 
occurs cease to cause sputtering of this atom, but the points out- 
side of this area still eject no. 19. At high energy the primary 
is reflected from the third layer atoms 59 and 64, and takes no. 19 
with it as it comes back out the front. 

At 10 keV atoms 13, 38 and 40, in addition to nos. 14 and 15, 
are sputtered often. Atom 13 is sputtered as a result of the primary 
driving 14 behind 13, whence the latter acquires motion out and to 
the left, rebounding from no. 12 before emerging. Atom 38 is sputter- 
ed by an identical mechanism, except that the primary, after inter- 
acting with no. 39 gets behind no. 38, and the latter interacts with 
37 before it comes out of the crystal. In an analogous manner no. 

39 sputters atom 40 . 

An interesting, but not surprising, phenomenon is the number 
of atoms sputtered from deep (third or greater) layers as a function 
of primary energy. There are none sputtered at 1 keV, one at 3 keV, 
four at 5 keV, and six at 10 keV. This behavior would be expected 
since the primary penetrates further with increased energy and 
transfers momentum to the deeper atoms. 
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As mentioned previously in section 2, the two potential approach 
resulted in a sputtering ratio curve of the proper shape. To adjust 
this curve to the ordinates obtained experimentally [24] , a "thresh- 
old energy" notion is employed. An arbitrary value is chosen, and 
the assumption made that any atom with kinetic energy less than this 
amount (in the direction normal to the surface) will not escape the 
crystal. A family of curves result (see fig. 5) from which an 
estimate can be made of the binding energy of atoms to a crystal 
surface. Thus 4 eV is the value estimated for the copper (110) 
surface since this figure yields a curve with the best fit to 
experimental data. 

An analysis of the angular distribution of the sputtered atoms 
leads to some extremely interesting results. At 1 keV primary 
energy one third of the ejected atoms leave the surface with equal 
velocities in the "x" and "z" directions. These velocities range 
from one-half to equal the "y" velocity. That is, a large percent- 
age of atoms emerge approximately in the ^111^ direction, but nearer 
to the surface normal. One-sixth of the atoms emerge in the ^110^ 
direction while the remainder are randomly oriented. Thus the model 
predicts a deposit spot pattern consisting of a central spot, four 
^111^ spots and considerable haze. 

At 10 keV one-third of the atoms emerge in the ^110^ direction, 
one-fifth in approximately the ^111^ direction, and one-fifth near 
the ^100^ direction. At this energy the model predicts the central 
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spot, four ^111^ spots, plus two ^100^ spots which may merge with 
the central spot, and less haze than obtained at 1 keV. 

2. The (111) Sample 

As an example of the sputtering mechanism, consider atom 19 , 
which is ejected most often. This atom has motion imparted by the 
primary "up’ 1 and towards the right (see fig. 3). The collision 
with atom 25 alters its momentum vector towards the front surface. 

The sputtering ratio curve is shown in figure 7. Although the 
shape of the curve approaches the experimental data for 6 eV threshold 
energy, the magnitudes disagree significantly. Inspection of the 
figure shows that different threshold energies do not yield a family 
of curves, but that the curve shape is altered with changes in 
energy. 

The deposit pattern for this orientation consists of a central 
spot with three large surrounding spots corresponding to the ^011^ 
directions . 



3. The (100) Sample 

Analysis of the paths of sputtered atoms reveals that the 
ejection mechanism is consistent with the other two orientations. 
Figure 9 shows that the atoms near the impact area are ejected most 
often. 

The sputtering ratio curve is shown in figure 6. There is close 
agreement with the experimental data in the 1-5 keV range. At 

10 keV the sputtering ratio is too high, indicating that the potential 
function may be too hard at higher energies. 
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The deposit pattern for this orientation consists of a heavy 
concentration of atoms in the center, with lesser concentrations 
in diagonally opposite quadrants. 
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4. Discussion 



A. The (110) Surface 

As mentioned in section 3, the majority of the atoms which are 
sputtered come from the region around the impact area* This obser- 
vation indicates a model which accounts for sputtering as the result 
of the direct transfer of momentum from the primary ion to the lattice 
atoms, rather than by chaining. The mechanism whereby crystal atoms 
are ejected supports this view. Although there are differences in 
this mechanism at high and low energies, there is no evidence in 
this study that momentum focusing plays a role in the sputtering 
process. 

The threshold energy necessary to achieve agreement with the 
experimental sputtering ratio curve is considered to be quite reason- 
able. The value of 4 eV is typical of the heat of sublimation, which 
is less than the energy required to displace an atom from its lattice 
site. The displacement energies which have been determined are 
applicable to interior atoms, while surface atoms are expected to 
be less tightly bound. 

The model successfully predicts the shape of the sputtering ratio 
curve. Discrepancies between experimental and theoretical magnitudes 
can be attributed to the approximations of the potential function 
in the theory or the modification of the crystal during irradiation 
experiments . 

The deposit spot pattern observed by Nelson, Thompson, Farmery 



16 



and Hall [p'Q is also successfully predicted. Although the experiment- 
al group was working at higher energies, deposit patterns above 10 
keV are expected to have the same general features. At 1 keV this 
study predicts a spot pattern similar to those obtained by Southern 
et $ al. [25]. This group worked with a ^011^ sample but crystall- 
ographic considerations show that the <110^ should have similar 
patterns. The model agrees qualitatively with the experimental 
results of the latter group in two important respects; (1) the non- 
central ejection spots are closer to the surface normal than would 
be expected by consideration of crystal geometry, and (2) the neigh- 
boring rows influence the motion of the ejected atoms. 

B. The (111) Surface 

The sputtering mechanism is similar to the (110) sample in that 
momentum focusing is not evident. The atoms are ejected as the 
result of sequences of first and second layer collisions. 

The estimated binding energy is taken to be 6 eV, the value 
which yields the proper shape of the sputtering ratio curve. Dis- 
crepancies in magnitude between the model and experiment are consider- 
ed to be the result of failure of the potential function parameters. 
"Hardening’* the potential for future studies of this orientation may 
give better agreement. 

Selection of different binding energies have the greatest effect 
on the magnitude of the sputtering ratio at lower energies. This 
indicates that lower primary energies result in ejection of lower 
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energy atoms. 

The deposit pattern experimentally observed by Southern et, al. 
L 2 d is successfully predicted. 

C. The (100) Surface 

Results of bombardment of this surface also indicate that 
momentum focusing is not evident. Sputtering is the result of first 
or second layer collisions. 

The estimated binding energy is 2.6 eV. Additional investigation 
in the 5-10 keV range is being conducted to bring the sputtering 
ratio curve into closer agreement with experimental observations. 

The shape of the curve is sensitive to the threshold energy, which 
has the greatest effect on magnitude at 3 keV (see fig. 6). The 
possibility exists, in nature, that the binding energy of the crystal 
is a function of the amount of disturbance, and hence a function of 
primary energy. 

The deposit pattern corresponds to "half" of the pattern 
observed by Southern et al. t 25 !- That is, there is a central spot 
and two ^011^ spots. Bombarding the crystallite in a mirror image 
impact triangle may yield the other two ^011^ spots. 
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5. Conclusions 



This computer simulation program is considered adequate in most 
respects. The model predicts the proper shape of the sputtering 
ratio curves, the angular distribution of ejected atoms, binding 
energies in the neighborhood of expected values, and magnitudes of 
the sputtering ratio for two of the three orientations. 

Momentum focusing does not enter into the sputtering mechanism 
for the argon-copper system in the 1- to 10-keV range. However, 
there are cases of chain transmission of momentum of surface atoms 
parallel to the surface. 

There is a wealth of experimental data available which requires 
theoretical exploration. Further work using this model should 
include study of the energy distribution of ejected atoms, sputtering 
ratio as a function of primary incoming direction, inclusion of 
thermalized motion of lattice atoms, and investigation of transmission 
sputtering, all for additional systems of elements and other energy 
ranges. The computer program should be expanded to include a graph 
plot output of the sputtering ratio curve and a plot of the angular 
and energy distributions of the ejected atoms. 

The computer programming assistance of LT Peter D.Stogisis 
gratefully acknowledged. I am indebted to Prof. Don E. Harrison Jr. 
for his guidance and assistance in all phases of the work. 

Note in support of argument . Recent work on argon sputtered 
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lithium by Stein, Malakhof and Smith QjfQ led the authors to 
suggest that, at low energy, there is small likelihood of extensive 
penetration of the lattice by the primary. Hence, there is little 
chance for a large number of n-body collisions yielding a focused 
momentum change to take place. The lobes of the distribution 
pattern of the sputtered atoms are believed to be created by a 
limited number of n-body collisions near the surface. (This 
effect should be even more pronounced for copper since the mass 
of the target atoms is considerably greater than lithium) . 
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FACE CENTERED CUBIC (100) 
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Face Centered Cubic (III) 
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APPENDIX I 



CRYSTAL GEOMETRY 

A. Lattice Dimmensions 



Orientation 


Dimensions (L.U.) 


X 


Y 


Z 


(no) 


9 x .707 = 6.36 


4 x .707 = 2.83 


8 x 1.0 - 8.0 


(100) 


6 x 1.0 - 6.0 


4 x 1.0 = 4.0 


8 x 1.0 = 8.0 


(HI) 


9 x .707 = 6.36 


4x1.155 - 4.62 


5.67x1.22= 6.94 



1 L.U. = 1.807 1 



B. Atomic Numbering System 





Y-plane number 


Orientation 


1 


2 


3 


4 


5 


(no) 


2-26 


27-46 


47-71 


72-91 


92-116 


(100) 


2-33 


34-64 


65-96 


97-127 


128-159 


(111) 


2-31 


32-61 


62-91 


91-121 


122-151 
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APPENDIX III 



PROGRAM DISCUSSION 



A. Brief Explanation of the Program. 

The descriptions below pertain to the sections of the program as 
outlined in appendix II. 



Section 1: 



Section 2s 
Section 3s 
Section 4s 

Section 5s 



Section 6s 



Allocates space in computer cells for arrays of various 
properties of the atoms; establishes input and output 
formats; establishes the Born-Mayer form for both 
potentials . 

Defines various constants and variables; reads input data. 
Constructs the crystallite. 

Commences the bombardment by assigning the primary its 
proper direction and velocity. 

Force double iteration procedure. This section calculates 
the forces on the atoms by use of a triangular matrix, 
then later recalculates the forces after the atoms have 
been displaced by the original forces. The average of 
these forces is then applied to the initial positions 
of the atoms to determine their displacements. 

Determines the tentative positions of the atoms under 
the forces calculated in section 5, then returns the 
program to section 5 for force recalculation. 
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Section 



Section 



Section 



Section 



Section 



Section 



Section 



7; Establishes positions of the atoms under the average force. 
Calculates velocities of all atoms. 

8: Tests the current time step to see if output is called 

for. If affirmative, prints output headings. 

9; Calculates atomic separations from the positions deter- 
mined in section 7. Calculates potential energies based 
upon these separations and the Born-Mayer parameters. 

10: Computes kinetic energies of the atoms from their velocities. 

Calculates the total binding energy of escaping atoms. 

11: If the time step is such that output is called for this 

section prints displacement, velocity, kinetic and 
potential energies of each atom with a total energy 
greater than 0.25 eV. 

12: Indexes any atoms which have escaped from the crystal so 

that they will be excluded from future "do” loops (an 
attempt to reduce computer running time). Recalculates 
the time step interval on the basis of the most energetic 
atom still in the crystal. Tests the total potential 
energy and time step interval against values selected 
to be indicative of the end of a run. Continues the run 
if tests are negative. 

13: Prints the status report at the end of the run, which 

includes the position, velocity, kinetic and potential 
energies of each atom. Returns the program to the input 
step for reception of the next impact point data. 
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B. Terms not defined in Leeds* Program ( 34 ] . 

PPOTF(X) . PFORCF(X) . PFXA . PFRC and PPTC are the argon-copper functions 
which correspond to the copper-copper functions P0TF(X),F0RCF(X), 
FXA,FRC and PTC. 

LCUT indexes any atoms which have been ejected from the crystal as 
explained under section 12, appendix II A. 

NLDTI is the number of time steps between recalculation of DT, the 
time step interval. 

NL is an index used to test if NLDTI time steps have elapsed since 
DT was last recalculated. 

QUIT is the value of total potential energy chosen to test for 
program shutdown. 

SLOW is the value of the time step interval chosen to test for 
program shutdown. 

BENERGY is the binding force assigned to the crystal surface, in eV. 
CVB is the conversion factor to make BENERGY consistent with units 
employed within the program. 

BRY is the cumulative distance from the front surface of all atoms 
that have been ejected. 

BF is the cumulative binding energy of these ejected atoms. The value 
is determined by a force x distance calculation (BENERGY x CVB x BRY). 
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APPENDIX IV 
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